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change in the number of zero modes as a function of chemical potential. 
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Quantum chromodynamics (QCD) at nonzero baryon density is important for the study of 
relativistic heavy-ion coUisions, neutron stars, and the early universe The effects of such a 
density are investigated by introducing a quark chemical potential pL in the QCD Dirac operator. 

For jU 7^ the Dirac operator loses its hermiticity properties and its spectrum moves into the 
complex plane. This causes a variety of problems, both analytically and numerically. Lattice sim- 
ulations are the main source of nonperturbative information about QCD, but at 7^ they cannot 
be performed by standard importance sampling methods because the measure of the functional 
integral, which includes the complex fermion determinant, is no longer positive definite. 

A better analytical understanding of QCD at very high baryon density has been obtained by 
a number of methods [||], and the QCD phase diagram has been studied in model calculations 
based on symmetries [Q]. Chiral random matrix theory (RMT) [|5|], which makes exact analytical 
predictions for the correlations of the small Dirac eigenvalues, has been extended to /X 7^ [^], 
and a mechanism was identified [0] by which the chiral condensate at /i 7^ is built up from the 
spectral density of the Dirac operator in the complex plane, in stark contrast to the Banks-Casher 
mechanism at /X = 0. This mechanism and its relation to the sign problem was also discussed by 
Splittorff in a plenary talk at this conference 

A first comparison of lattice data with RMT predictions at 7^ was made in Ref. using 
staggered fermions. One issue with staggered fermions is that the topology of the gauge field is 
only visible in the Dirac spectrum if the lattice spacing is small and various improvement and/or 
smearing schemes are applied [[l0[]. To avoid these issues, we would Uke to work with a Dirac op- 
erator that implements a lattice version of chiral symmetry and has exact zero modes at finite lattice 



spacing [11]. The overlap operator [12] satisfies these requirements at /i = 0. In the following, we 
show how the overlap operator can be modified to include a nonzero quark chemical potential' [jl|]. 
We then study the spectral properties of this operator as a function of /i and compare data from 
lattice simulations with RMT predictions. As we shall see, the overlap operator has exact zero 
modes also at nonzero pL, which allows us, for the first time, to test predictions of non-Hermitian 
RMT for nontrivial topology. 

We begin with the well-known definition of the Wilson Dirac operator including a chemi- 



cal potential /i [15], 
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Dwipi) = 1 - (7;.+ + 7;.-) - K (e^'T^ + e-^'T^-) , 

(r±),, = (1 ± 7v)t/±v W5j.,±v , (1) 

where k = 1/ {2mw + 8) with the Wilson mass m^, the U G SU(3) are the lattice gauge fields, and 
the 7v are the usual Euclidean Dirac matrices. Unless displayed explicitly, the lattice spacing a is 
set to unity. 



The overlap operator is defined at /i = by [ 12] 



Dov(O) = l + 75£(75£>w(0)), (2) 



' See also Ref. [ |l3| | for a perfect lattice action at /I 7^ and Ref. [ |l4| ] for an overlap-type operator at /i 7^ in 
momentum space. 
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where e is the matrix sign function and 75 = 71727374. m^/ must be in the range (—2,0) for Dov(O) 
to describe a single Dirac fermion in the continuum. The properties of Dov(O) have been studied 
in great detail in the past years. In particular, its eigenvalues are on a circle in the complex plane 
with center at (1,0) and radius 1, its nonreal eigenvalues come in complex conjugate pairs, and it 
can have exact zero modes without fine-tuning. Dov(O) satisfies a Ginsparg -Wilson relation [|l^] of 
the form 

{D,75} = D75D. (3) 

We now extend the definition of the overlap operator to 7^ 0. The operator D-w (0) in Eq. 
is 75-Hermitian, i.e., 75£>w(0)75 =d]^{0), and therefore the operator y^D-v/{0) in the matrix sign 
function is Hermitian. However, for /i 7^ 0, Dw{^) is no longer 75-Hermitian. Defining the overlap 
operator at nonzero /i by 

Dov(At) = 1 + 75e(750w(iU)) , (4) 

we now need the sign function of a non-Hermitian matrix. In general, a function / of a non- 
Hermitian matrix A can be defined by the contour integral 

f{A) = ^.j^dzf{z){zl-A)-\ (5) 

where the spectrum of A is enclosed by the contour T and the matrix integral is defined on an 
element-by-element basis. A more convenient expression can be obtained if A is diagonalizable. In 
this case we can write A = UAU^\ where U G Gl{N,C) with A'^ = dim (A) and A = diag(Ai, . . .,Xn) 
with A,- € C. Then, applying Cauchy's theorem to Eq. (|), /(A) = Uf{A)U-\ where /(A) is a 



diagonal matrix with elements /(A,). In particular, the matrix sign function can be defined by [|T7|] 

£(A) =[/sign(ReA)[/"^ (6) 

This definition ensures that £^(A) = 1 and gives the correct result if A is real. An equivalent 
definition is e(A) =A{A^y^/^ Eqs. (|) and (|) constitute our definition of Dov(At). The 

sign function is ill-defined if one of the A, lies on the imaginary axis. Also, it could happen that 
75Z)iv(/x) is not diagonalizable (one would then resort to a Jordan block decomposition). Both of 
these cases are only realized if one or more parameters are fine-tuned, and are unlikely to occur in 
realistic lattice simulations. 

It is relatively straightforward to derive the following properties of Dov(jtt): 

• It is no longer 75-Hermitian. 

• It still satisfies the Ginsparg- Wilson relation (|3|) because of £^(A) = 1. Thus, we still have a 
lattice version of chiral symmetry, and the operator has exact zero modes without fine-tuning. 

• Its eigenvalues not equal to or 2 no longer come in complex conjugate pairs, but every such 
eigenvalue A (with eigenvector comes with a partner A /(A — 1) (with eigenvector 751^). 

• Its eigenvectors corresponding to eigenvalues or 2 can be arranged to have definite chiraUty. 

We now turn to our (quenched) lattice simulations. The computation of the sign function of 
a non-Hermitian matrix is very demanding. We are currently investigating various approximation 
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Figure 1: Spectrum of DovifJ-) for jj. = and jj. = 0.3 for a typical configuration. The figure on the right is 
a magnification of the region near zero. 



schemes, but in this initial study we decided to compute the sign function and to diagonalize Dovip.) 
exactly using LAPACK. For the comparison with RMT we need high statistics, which restricts us 



to a very small lattice size. We have chosen the same parameter set as in Ref. [19] to be able to 
compare with previous results at = 0. The lattice size is V = 4^, the coupling in the standard 
Wilson action is j8 = 5. 1, the Wilson mass is mw = —2, and the quark mass is = 0. 

In Fig. |l| we show the spectrum of Dov(iU) for a typical configuration for /i = and n = 0.3. 
As expected, we see that the eigenvalues move away from the circle as n is turned on. Another 
observation is that the number of zero modes of Dov(jU) for a given configuration can change as a 
function of jJ., see Figure]^ This can be understood from the relation between the anomaly and the 



index of Dqv |2C, 21|, 



-tr(75Z)ov) =2index(Do 



(7) 



which we can show to remain valid at |i 7^ 0. Using tr(75Dov) = tr[e{Y5Dw)] and the fact that 
the eigenvalues of the sign function are +1 or —1, one has index(Dov) = — n^)/2, where 
denotes the number of eigenvalues of y^DwilJ') with real part ^ 0. Therefore the number of zero 
modes for a given configuration is determined by the difference of the number of eigenvalues of 
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Figure 2: Distribution of the number of zero mores of Dov(/i) for jj, = 0,0.1,0.2,0.3, 1.0, and number of 
configurations (right table). 
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Figure 3: Scatter plot of the projected eigenvalues A of Dov(/^) for /i = 0.1 (left pane) and pL — 0.3 (right 
pane). The eigenvalues are projected using A — > 2A/(2 — A). At /i = this mapping projects the eigenvalues 
from the GW-circle to the imaginary axis. At /i 0, the same mapping projects the eigenvalues onto a band 
parallel to the imaginary axis. As /i is increased the eigenvalues spread further into the complex plane. 



y^Dw lying left and right of the imaginary axis. As pi changes, an eigenvalue of J^Dw can move 
across the imaginary axis. As a result, index(Dov) changes by 1, which explains the observation. 
We believe that this is a lattice artefact which will disappear in the continuum limit. 

The spectral density of Dov(jU) is given by p°^(Ar,A;) = (L/t^C'^ — ^k)) with X = Xr + iXi, 
where the average is over configurations. The distribution of (projected) eigenvalues in the complex 
plane is shown in Figure |3[ The claim is that the distribution of the small eigenvalues of Dov{jx) is 
universal and given by RMT The chiral RMT model for the Dirac operator is ^ 




where W is a complex matrix of dimension « x {n + v) with no further symmetries (we take v > 
without loss of generality). The matrix in Eq. (^ has v eigenvalues equal to zero. The spectral 
correlations of Drmt(m) on the scale of the mean level spacing were computed in Refs. [22, 23, 2^]. 
In the quenched approximation, the result for the microscopic spectral density reads 



where z = x + iy = AZV, / and K are modified Bessel functions, and a = }j}f^V. Z and /;;: are 
low-energy constants that can be obtained from a two-parameter fit of the lattice data to the RMT 
prediction, Eq. (^). (The normalization is fixed by / dxdy p*^^ {x,y) = \2V and does not introduce 
another parameter.) For x <C a, Eq. (Rl) becomes radially symmetric 



pT^''{x,y)^:^KA^)lv{^) (10) 

with 1^ = |zp/4a, and the fit only involves the single parameter ^/ fn- 

In Fig. ^ we compare our lattice data to the RMT prediction. We display various cuts of the 
eigenvalue density in the complex plane as explained in the figure captions. The data agree with 
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Figure 4: Density of the small eigenvalues z = x + /y = A EV of Dov(/x) in the complex plane (after projecting 
A ^ 2A/(2 - A)) for (from left to right) pi = 0.1,0.2,0.3, 1.0. The histograms are lattice data for v =0 and 
V = 1, and the solid lines are the corresponding RMT prediction of Eq. (^, integrated over the bin size. 
Top: cut along the imaginary axis, middle: cut along the real axis, bottom: cut parallel to the real axis at 
yc- Vertical lines indicate the fit interval. For ^ = 1.0 the distribution of the small eigenvalues is radially 
symmetric up to |^ | 0.7, with ^ = |z|^/4a, and therefore only the ratio T./ fn can be determined from a fit 
to Eq. (p^. In the rightmost bottom plot the data are integrated over the phase. 
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Table 1: Fit results for E and fj^. For ii — \ .Q only E//;r 
can be determined (see Figure 

the RMT predictions within our statistics. £ and fji were obtained by a combined fit to the v = 
and V = 1 data for all three cuts and are displayed in Table [I[ (These numbers have no physical 
significance at j8 = 5.1.) 

In summary, we have shown how to include a quark chemical potential in the overlap operator. 
The operator still satisfies a Ginsparg- Wilson relation and has exact zero modes. The distribution 
of its small eigenvalues agrees with predictions of non-Hermitian RMT for trivial and nontrivial 
topology. Our initial lattice study should be extended to weaker coupling, larger lattices, and 
better statistics. Work on approximation methods to enable such studies is in progress. For small 
volumes, reweighting with the fermion determinant should allow us to test RMT predictions for 
the unquenched theory 
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This work was supported in part by DFG. The simulations were performed on a QCDOC 



machine in Regensburg using USQCD software and Chroma [|27|], and GotoBLAS optimized for 
QCDOC. 
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